SYSTEM AND METHOD TO INTERPRET DISTRIBUTED 
TEMPERATURE SENSOR DATA AND TO DETERMINE A FLOW 
RATE IN A WELL 

5 DESCRIPTION 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[Para 1] This application claims the benefit under 35 USC 
1 1 9(e) of US Provisional Patent Application Serial Number 
1 0 60/5331 88 filed 30 December 2003 and US Provisional Patent 
Application Serial Number 60/536059 filed 1 3 January 2004. 

BACKGROUND 

[Para 2] The invention generally relates to a system and 
1 5 method to interpret distributed temperature sensor data. More 
specifically, the invention relates, first, to a system and method 
to automatically interpret distributed temperature sensor data 
from a specific wellbore event, such as the thermal events 
typically present in a gas lift well, second, to a system and 
20 method to determine flow rates from temperature data 

obtained from a gas lift wellbore, and third, to a system and 
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method to optimize production flow rates from temperature 
data obtained from a gas lift wellbore. 

[Para 3] Distributed temperature sensors, such as Sensor 
Highway Limited's DTS line of fiber optic distributed 
5 temperature sensors, have been used to measure the 

temperature profile of subterranean wellbores. In the DTS 
systems, an optical fiber is deployed in the wellbore and is 
connected to an opto-electronic unit that transmits optical 
pulses into the optical fiber and receives returned signals back 

1 0 from the optical fiber. Depending on the type of wellbore and 
on the service or completion, the optical fiber may be deployed 
in a variety of ways, such as part of an intervention service, 
permanently inside of a tubing (such as a production tubing), 
or permanently installed in the annulus between the borehole 

1 5 wall and the tubing. The signal reflected from the optical fiber 
and received by the opto-electronic unit differs depending on 
the temperature at the originating point of the reflected signal. 

[Para 4] Sensor Highway's DTS system utilizes a technique 
called optical time domain reflectometry ("OTDR"), which 
20 detects Raman scattering to measure the temperature profile 
along the optical fiber as described in U.S. Pat. Nos. 4,823,166 
and 5,592,282 issued to Hartog, both of which are 
incorporated herein by reference. For purposes of 
completeness, OTDR will now be described, although it is 
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understood that OTDR is not the only way to obtain a 
distributed temperature measurement (and this patent is 
therefore not limited to OTDR). 

[Para 5] In OTDR, a pulse of optical energy is launched into 
5 an optical fiber and the backscattered optical energy returning 
from the fiber is observed as a function of time, which is 
proportional to distance along the fiber from which the 
backscattered light is received. This backscattered light 
includes the Rayleigh, Brillouin, and Raman spectra. The 

1 0 Raman spectrum is the most temperature sensitive with the 
intensity of the spectrum varying with temperature, although 
Brillouin scattering and in certain cases Rayleigh scattering are 
also temperature sensitive. Generally, in one embodiment, 
pulses of light at a fixed wavelength are transmitted from a 

1 5 light source down the fiber optic line. Light is back-scattered 
along the length of the optical fiber and returns to the 
instrument. Knowing the speed of light and the moment of 
arrival of the return signal enables its point of origin along the 
fiber line to be determined. Temperature stimulates the 

20 energy levels of molecules of the silica and of other index- 
modifying additives - such as germania - present in the fiber 
line. The back-scattered light contains upshifted and 
downshifted wavebands (such as the Stokes Raman and Anti- 
Stokes Raman portions of the back-scattered spectrum) which 
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can be analyzed to determine the temperature at origin. In 
this way the temperature along the fiber line can be calculated 
by the instrument, providing a complete temperature profile 
along the length of the fiber line. Different temperature 
5 profiles can also be obtained in time, thereby providing a time 
lapsed temperature profile along the entire length of the 
optical fiber. 

[Para 6] The temperature profiles that are obtained from 
distributed temperature sensors such as the DTS can then be 

1 0 used by operators to, among others, measure flow rate, 
identify the presence and location of leaks, or identify the 
extent and success of an injection operation. However, the 
temperature profiles obtained from distributed temperature 
sensors such as the DTS generate a very large amount of data 

1 5 per time profile. This data is currently typically reviewed 
manually, at least at some point during the analysis. 
Reviewing this data manually in order to analyze and extract 
value from the data is a time consuming and highly specialized 
operation. 

20 [Para 7] For instance, the temperature profiles generated 

from distributed temperature sensors are very useful in gas-lift 
operations. Gas expands abruptly where it enters production 
tubing. This expansion produces significant cooling through 
the Joule-Thomson effect. Consequently, the temperature 
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profiles can reveal where, when, and to what extent gas is 
injected (i.e. the location and operation of a gas lift valve). 
However, the temperature profile often fluctuates in gas lift 
wells. Although the temperature change at an injection point 
5 may be several degrees Centigrade, the presence of 

fluctuations, the exceedingly high number of temperature data 
points, and the broad temperature trend in the well may 
obscure the change. Thus, it takes an operator a substantial 
amount of time to manually identify the sections of the 

1 0 temperature profile that contain valuable information and to 
then remove or suppress the background or non-relevant 
temperature phenomena from the valuable information in the 
temperature profile. Manual analysis introduces subjectivity, 
cannot be automatically integrated with use of other 

1 5 algorithms, and may provide inaccurate analysis due to noise 
or temperature trends that obscure the signal to a human 
operator. Furthermore, the ability to obtain production flow 
related information from the distributed temperature sensor 
data or other temperature data has been limited. 

20 [Para 8] Thus, there is a continuing need to address one or 
more of the problems stated above. 
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[Para 9] A distributed temperature sensor is deployed in a 
wellbore and is functionally connected to a processor. The 
processor receives the data from the distributed temperature 
sensor and automatically processes the data to highlight 
5 valuable information to the user relating to the relevant well, 
completion, or reservoir. In one embodiment, the distributed 
temperature sensor and processor are utilized in a well 
containing a gas-lift system, wherein the processor highlights 
valuable information to the user pertaining to the gas-lift 

1 0 system. A well model enabled by the processor enables the 

determination of a produced fluid flow rate in the well having a 
gas lift system. Temperatures are measured within the well to 
obtain a temperature profile, and this profile is processed 
according to the well model. The well model relates thermal 

1 5 characteristics, e.g. thermal decay and/or amplitude of a 
thermal discontinuity at an injection point, to flow rate. 

[Para 1 0] Advantages and other features of the invention will 
become apparent from the following drawing, description and 
claims. 

20 

BRIEF DESCRIPTION OF THE DRAWING 
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[Para 1 1 ] Fig. 1 is a schematic of a well completion utilizing 
the present Invention, including a distributed temperature 
sensor and a processor. 

[Para 1 2] Fig. 2 is a schematic of one embodiment of the 
5 processor. 

[Para 1 3] Fig. 3 is a flow chart of one algorithm that may be 
performed by the processor. 

[Para 1 4] Fig. 4 is a schematic of a well completion with a gas 
lift system utilizing the present invention. 

1 0 [Para 1 5] Figs. 5-1 0 are flow charts of different parts or 
embodiments of algorithms that may be performed by the 
processor. 

[Para 1 6] Figs. 1 1 -1 6 are plots that may be presented to an 
operator to illustrate the results of the operations performed 
1 5 by the processor. 

[Para 1 7] Figs. 1 7-22 are flow charts representing algorithms 
and specific portions of those algorithms according to 
alternate embodiments of the present invention. 

[Para 1 8] Figs. 23-24 are flow charts representing an 
20 algorithm for optimizing gas-lift performance of a well. 

DETAILED DESCRIPTION 
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[Para 1 9] Figure 1 illustrates an embodiment of a system 1 0 
that is the subject of this invention. A wellbore 12 extends 
from the surface 1 1 into the earth and intersects a formation 
14 that contains fluids, such as hydrocarbons. The wellbore 
5 1 2 may be cased. A tubing 1 6, such as a production tubing, 
extends within the wellbore 1 2. A packer 1 5 provides a seal 
and isolates the formation 14 from the region above the 
packer 1 5. Depending on whether the wellbore 1 2 is used as 
an injector well or as a producing well, fluid is either injected 

1 0 into the tubing 1 6 and into the formation 1 4 or fluid is 

produced from the formation 1 4 and into the tubing 1 6. In 
either case, fluid enters or exits the tubing 1 6 through flow 
paths in the tubing 1 6, such as the ports 1 8 illustrated in 
Figure 1 . The injection and production of fluids may also be 

1 5 aided by artificial lift mechanisms, such as pumps or gas lift 
valves. Perforations (not shown) may also be made in the 
wellbore 1 2 at the formation 1 4 in order to facilitate the flow 
of fluids into or out of the formation 14. 

[Para 20] The system 1 0 includes a distributed temperature 
20 sensor system 20 and a processor 22. The sensor system 20 
comprises an optical fiber 24 deployed in the wellbore 1 2 and 
an opto-electronic unit 26 typically but not necessarily located 
at the surface 1 1 . The optical fiber 24 is connected to the unit 
26 and can be used to measure temperature simultaneously at 
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multiple depths. In one embodiment, the optical fiber is 
deployed within a control conduit 28, such as a 0.25" control 
line. The conduit 28 may be attached to the tubing 1 6. 
Although the conduit 28 is shown attached to the exterior of 
5 tubing 1 6, conduit 28 (and optical fiber 24) may instead be 
inside of tubing 16 or may be cemented in place to the outside 
of the casing (not shown). In one embodiment, the optical 
fiber 24 is injected into the conduit 28, which may also be u- 
shaped, byway of fluid drag, as disclosed in U.S. Pat. No. Re 
10 37283, which patent is incorporated herein by reference. The 
optical fiber 24 also may be implemented as a temporary 
distributed temperature sensor installation or as a slickline 
distributed temperature sensor system. 

[Para 21] As previously disclosed, the unit 26 launches optical 
1 5 pulses into the optical fiber 24 and backscattered light is 
returned from the optical fiber 24. The backscattered light 
signals include information which can provide a temperature 
profile along the length of the optical fiber 24. For the 
configuration of Figure 1 , the temperature profiles generated 
20 by the system 1 0 may be used to detect whether fluids are 
flowing from the formation 1 4 and into the tubing 1 6 or to 
detect the extent and success of an injection operation from 
the tubing 1 6 and into the formation 1 4. 
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[Para 22] Processor 22 automatically analyzes the temperature 
profile data to minimize or remove any non-relevant 
temperature "noise" and to focus on tlie data points or sections 
tliat contain valuable information. As will be described, the 
5 processor 22 may also be programmed by an operator to " 
identify" particular temperature signals that typically 
correspond to a particular downhole event having an inflow of 
cooler fluid, e.g. gas, into a flowing stream, e.g. a flowing 
stream of oil or oil, water and gas mixtures. These types of 

1 0 events indicate, for example, the location of a gas lift valve, a 
hole in the production tubing, general wellbore completion 
tool leaks (e.g. packer leaks, sliding sleeve leaks, collar leaks) 
or the inflow of fluids from a formation that are cooler than the 
fluid flowing in the wellbore. The cooler temperatures typically 

1 5 are due to Joule-Thompson expansion of the inflowing fluid at 
or near the inflow point and indicate the magnitude of the 
inflow to the continuous flow stream and whether it is a 
continuous or transient event. 

[Para 23] The processor 22 is connected to the unit 26 by way 
20 of a communication link 30. The communication link 30 can 
take various forms, including a hardline, e.g. a direct hard-line 
connection at the well site, a wireless link, e.g. a satellite 
connection, a radio connection, a connection through a main 
central router, a modem connection, a web-based or internet 
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connection, a temporary connection, and/or a connection to a 
remote location such as the offices of an operator. The 
communication linic 30 may enable real time transmission of 
data or may enable time-lapsed transmission of data. The 
5 data transmission and processing allow a user to monitor the 
wellbore 1 2 in real time and take immediate corrective action 
based on the data received or analysis performed. In other 
words, processor 22 is able to process the data as it is 
received, enabling a controller/operator to make real-time 
10 decisions. 

[Para 24] Processor 22 may be a portable computer that can 
be removably attached from the unit 26. With the use of a 
portable computer, a user may analyze various wellbores while 
using a single computer system. Processor 22 may be a 
1 5 personal computer or other computer. 

[Para 25] Figure 2 illustrates in block diagram form an 
embodiment of hardware that may be used as the processor 
22 and to operate the representative embodiment of the 
present invention. The processor 22 comprises a central 
20 processing unit ("CPU") 32 coupled to a memory 34, an input 
device 36 (i.e., a user interface unit), and an output device 38 
(i.e., a visual interface unit). The input device 36 may be a 
keyboard, mouse, voice recognition unit, or any other device 
capable of receiving instructions. It is through the input device 
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36 that the user may make a selection or request as stipulated 
herein. The output device 38 may be a device that is capable 
of displaying or presenting data and/or diagrams to a user, 
such as a monitor. The memory 34 may be a primary memory, 
5 such as RAM, a secondary memory, such as a disk drive, a 
combination of those, as well as other types of memory. Note 
that the present invention may be implemented in a computer 
network, using the Internet, or other methods of 
interconnecting computers. Therefore, the memory 34 may be 

1 0 an independent memory 34 accessed by the network, or a 
memory 34 associated with one or more of the computers. 
Likewise, the input device 36 and output device 38 may be 
associated with any one or more of the computers of the 
network. Similarly, the system may utilize the capabilities of 

1 5 any one or more of the computers and a central network 
controller. Therefore, a reference to the components of the 
system herein may utilize individual components in a network 
of devices. Other types of computer systems also may be 
used. Therefore, when reference is made to "the CPU," "the 

20 memory," "the input device," and "the output device," the 

relevant device could be any one in the system of computers or 
network. 

[Para 26] Figure 3 shows in flow chart form the operations 
performed by the processor 22. In the first step 40, the 
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processor 22 obtains the distributed temperature sensor data 
(i.e. tlie temperature profiles) from the unit 26 via the 
communication link 30 as previously disclosed. In step 42, the 
processor 22 processes the data to remove noise and/or focus 
5 on significant events to thereby extract the valuable 
information from the data. In step 44, the valuable 
information is provided as an output In the format chosen by 
the user through the output device 38. 

[Para 27] The process data step 42 can tal<e on a variety of 
1 0 forms, depending on the desire of the operator and on the 

configuration of the wellbore being analyzed (i.e. gas lift, water 
injection, producer, horizontal). In one embodiment, the 
process data step comprises the use of an algorithm to 
process the temperature profile data to remove noise from the 
1 5 data and/or focus on significant events. Generally, the 

algorithms that may be used to achieve these functions include 
the removal of low order spatial trends (e.g. a polynomial In 
depth can be fit to each temperature profile and the resulting 
function can be subtracted from the profile), a high-pass filter 
20 (such as one that removes low spatial frequencies like a sixth- 
order, zero-phase Butterworth filter), the differentiation of 
data with respect to an Independent variable (such as depth), 
low pass filters, matched filters (functions with shapes similar 
to what Is expected In the data), adaptive filters, wavelets. 
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background subtraction, Bayesian analysis, and model fitting. 
These algorithms can be applied to the data individually, or in 
combination. For example, filtering can identif/ important 
regions of the data and then trend removal can be used for 
5 further processing. Moreover, the algorithms can be applied in 
measured depth or in time. It should be noted the algorithm 
may be applied to other applications, such as detection of 
carbon dioxide or steam flood in production wells and to 
identify other events having a large Joule-Thompson effect. 

1 0 [Para 28] An example of how the model fitting algorithm may 
be used to analyze a gas-lift well will now be described with 
reference to Figure 4. Figure 4 illustrates the wellbore 1 2 
including a gas lift system 50 disposed therein. Gas lift system 
50 may comprise at least one gas lift valve 52 disposed on 

1 5 tubing, such as production tubing 1 6. Gas lift is a common 
method for providing artificial lift to oil wells and involves 
injecting gas 54 from a source 56 into the annulus 58 through 
the valves 52 and into tubing 1 6. The valves 52 are typically 
pressure-controlled, with only the deepest valve 52 open 

20 during normal operation. Shallower valves are opened to start 
the well flowing. The gas reduces the average density in the 
production column by displacing oil and water. Thus, the gas 
injection increases production by reducing the pressure at the 
bottom of the well. 
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[Para 29] The design of the gas-lift system 50 is matched to 
the productive capacity of the well. Design parameters include 
gas-Injection pressure and rate, tubing diameter, valve depths 
and operating pressures, and orifice diameters of the valves. 
5 However, equipment failures, changes In a well's In-flow 

capacity, or changes in water-cut can reduce the effectiveness 
of the gas-lift system. Because gas injection often causes 
large fluctuations In production, traditional production logging 
tools, which measure at each depth at a different time, can 
10 provide ambiguous data. Consequently, diagnosing problems 
in gas-lift wells is difficult. The time-lapsed temperature 
profiles generated by distributed temperature sensors are 
particularly suited and beneficial for this diagnosis. 

[Para 30] The use of the model fitting algorithm to analyze a 
1 5 gas-lift system 50 achieves the following: [1 ] it removes the 
Irrelevant aspects of the temperature profile data and 
suppresses noise, [2] It tolerates the rapid temperature 
fluctuations In space and time that are typical In gas-lift wells, 
[3] it tolerates the possibility that the gas signature may be 
20 limited to a small region or spread out over a large one, [4] it 
minimizes input from an operator thereby reducing training 
requirements and the staff time that must be devoted to 
processing, and [5] It processes the data rapidly making it 
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useful in temporary (and not only permanent) distributed 
temperature sensor installations. 

[Para 31 ] With the model fitting algorithm, a model of at least 
part of the wellbore or its performance is fit to the temperature 
5 profile data by adjusting parameters in the model. As 
illustrated in Figure 5, the process data step 42 in this 
embodiment comprises generating the model at step 46 and 
then fitting the model to the data at step 48. 

[Para 32] One embodiment of the fitting the model to the data 
1 0 step 48 is shown in Figure 6. In this embodiment, the model's 
function is calculated with initial parameter values that may be 
estimated at step 60. Then, at step 62, the parameter values 
are changed (unless it is the first iteration including the initial 
parameter values). At step 64, the model's function is 
1 5 calculated with the new parameter values. And, at step 66, the 
current model calculation and the subsequent model 
calculation are compared and the one that provides the better 
fit to the temperature profile data is stored. The iteration then 
continues until further modifications of model parameters no 
20 longer make the fit significantly better, at which point the 
initial model of that iteration is deemed to be the best fit for 
the temperature profile data. The definition of best fit may be 
preprogrammed by a user. 
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[Para 33] In one embodiment, the model comprises a 
comprehensive model for the physical gas-lift system. In 
another embodiment, the model comprises a 
phenomenological model. 

5 [Para 34] With respect to the use of a phenomenological 
model, the basic assumption of the model may be that gas 
injection causes a local perturbation of temperature and that 
the perturbation decreases exponentially in either direction 
from the point of injection. Because flowing fluids convect 

1 0 heat, the decay length is greater in the downstream (up the 
well) direction. The model is fit to a range around each valve. 
The model has four primary parameters for each valve: the 
depth of injection in the wellbore (i.e. approximate valve 
depth), the amplitude of the temperature effect (i.e. how much 

1 5 temperature difference is caused by the injection), and the 
decay length in each direction from the injection depth. The 
model also includes two secondary parameters for each valve, 
a slope and an intercept, to account for a linear background 
temperature variation. The model adjusts the parameters to 

20 match the data at each valve in each temperature profile. 

When the distance between valves is sufficient, the valves are 
treated independently because the effect at one valve caused 
by another is smooth and can be considered part of the 
background. Otherwise, valves that are close together may be 
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grouped for simultaneous analysis. The algorithm uses the 
Levenberg-Marquardt method, discussed in DW Marquardt, J. 
Soc. Industrial and Applied Mathematics, vol. 1 1 , p. 1 31 
(1 963), to solve the non-linear fitting problem. The algorithm 
5 also tests the fit at each valve in each profile for statistical 
significance. If a fit is not considered significant, the 
temperature amplitude is set to zero and other parameters are 
set to default values. 

[Para 35] A function that may be utilized for the 
1 0 phenomenological model described above is the following 
modification of one derived by Ramey (H.J. Ramey, "Wellbore 
heat transmission," J. Petroleum Technology, p427 (1 962)) for 
the thermal signature of fluids pumped down a well: 

Equation 1 

15 r,. =Aexp-l-y 

V 

wherein A is the amplitude of the temperature effect, zj\s the 
depth measured from the surface (a position depth variable), d 
is the approximate valve depth , a is the background intercept, 
and b is the background slope. If z/exceeds d, then h is the 
20 upstream (down the well) decay length and the downstream 

(up the well) decay length is ignored. If Z/is less than d, then // 
is the downstream (up the well) decay length and the upstream 
(down the well) decay length is ignored. 
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[Para 36] Figure 7 illustrates in flow chart form the algorithm 
used to solve the phenomenological model with use of 
Equation 1 . In the first step 70, the user is prompted for and 
the user inputs the approximate depth of each gas lift valve. 
5 The temperature profile data from the distributed temperature 
sensor is then read in step 72. in the third step 74, the 
algorithm selects a region for fitting for each valve. In this 
step 74, the algorithm sets up a window of data to analyze 
above and below the approximate valve depth inputted in step 

1 0 70 for each valve. In the next step 76, the statistical noise 
level in each region selected in step 74 is estimated for each 
time profile of temperature data. The model is then fit to the 
actual data for each valve and for each time profile of 
temperature data in step 78. Next, in step 80, the result of 

1 5 each fit is tested for statistical significance to ensure the fit is 
an actual and not a computer-created event. Lastly, in step 
82, the results are displayed in ways that are beneficial and 
valuable to the user. 

[Para 37] Although the algorithm illustrated in Figure 7 
20 prompts the user for the approximate valve depth in step 70, 
in another embodiment the algorithm automatically locates the 
location of each valve by interpreting and analyzing the data 
from the distributed temperature sensor. Given the typical 
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gas-lift valve signature, a matched filter may be used by the 
processor 22 to locate each valve. 

[Para 38] Step 76 (Estimate Statistical Noise Level For Each 
Region And Time Profile) is further illustrated in Figure 8. 
5 First, in step 84, regions very near to each valve are omitted 
and the remaining data are sub-divided by depth into groups. 
Next, in step 86, the linear trends are removed from the data 
in each of those groups. In step 88, the power spectrum of the 
relevant data is estimated in each group. And then, in step 90, 
1 0 the baseline in the power spectrum is estimated as the 

statistical noise level for that group. In step 91 , a smooth 
curve is fit to the noise level versus depth, and the noise level 
for the depth of each valve is taken from the smooth curve. 

[Para 39] Step 78 (Fit Model To Data For Each Valve In Each 
1 5 Time Profile) is further illustrated in Figure 9. In the first step 
92, the linear background coefficients from Equation 1 (a and 
b) are estimated, such as by selecting random starting points 
for each or by selecting a line with a given slope as a starting 
point. In the next step 94, the amplitude of the temperature 
20 perturbation caused by the injection is estimated by analyzing 
the temperature data. Then, in step 96, the six parameters of 
Equation 1 are adjusted to provide the best fit to the actual 
temperature data, such as by using the sum of squares of 
deviations method. 
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[Para 40] Step 96 (Adjust Parameters To Minimize Sum Of 
Squares Of Deviations) is further illustrated in Figure 1 0. in 
the first step 98, the model is computed using the initial values 
for the various parameters. Next, in step 1 00, constraints for 
5 valve depth and decay lengths are taken into account to ensure 
that such parameters are not iterated to be outside of certain 
ranges. For instance, in one embodiment, a constraint is 
placed on the valve depth parameter to ensure that it remains 
within the region selected in step 74 of Figure 7. And, a 

1 0 constraint is placed on the decay lengths to ensure that such 
values are always positive (not negative or 0). In the next step 
1 02, the values of each of the parameters are iterated with the 
goal of minimizing deviations. In step 1 04, the initial 
parameter estimates or values are perturbed or changed again, 

1 5 and preferably twice more, and the iteration process is rerun 
for each perturbation. This step 104 ensures that the global 
and not just a local minimum is generated as a result of the 
iteration process. In the last step 1 06, the best fit from each 
of the iteration sequences is selected as the best overall fit and 

20 the tolerance is reduced to make the final fit. 

[Para 41 ] It should be noted that the goal of step 80 (see 
Figure 7) is to reduce the number of non-events that are 
incorrectly identified as events ("false positives") and to reduce 
the number of events that are incorrectly ignored ("false 
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negatives"). In one embodiment, step 80 comprises comparing 
competing models. An appropriate competing model to the 
phenomenological model previously described is a low-order 
polynomial. Although the downstream temperature decay 
5 length for gas injection may range from a few meters to 
hundreds of meters, the upstream decay length should be 
limited to a few meters. Thus, an event has some temperature 
variation that occurs over a short distance. Non-events that 
the model may fit result primarily from the convection of 

1 0 temperature disturbances up the well. The sharp features of a 
temperature fluctuation smooth out as it travels. Since a low- 
order polynomial is smoother than the target model, it fits 
most convected features better. The residual variance, that is, 
the sum of squares of the differences between data and a 

1 5 model, is a common measure of the quality of a fit - the 

smaller the variance is, the better the fit. If the variance for the 
phenomenological model is smaller than that for the 
polynomial, the phenomenological model may be regarded as 
significant. In one embodiment, one may require that the 

20 variance of the polynomial model be larger by some fraction of 
the noise level. In another embodiment, the required fraction 
may be reduced when adjacent profiles have significant fits. 

[Para 42] In step 82 of Figure 7, the results are displayed to 
the user. Figure 1 1 illustrates one plot that may be useful to a 
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user. The plot is depth versus temperature and the points 108 
on the plot are the raw temperature profile data points 
obtained from the distributed temperature sensor near one of 
the gas lift valves at one particular time. The curve 1 1 0 on the 
5 plot is the curve derived from the phenomenological model 
and algorithm previous described that best fits the points 1 08. 
A review of this plot would enable a user to visualize the 
accuracy of the fitting, which in the case illustrated is good. 

[Para 43] The user may also want to view the pure 
1 0 temperature perturbation created by the injection at a specific 
valve without the background linear trend. In order to plot this 
pure perturbation, the background parameters (a and b) are 
removed from Equation 1 , giving: 

Equation 2 



and Eq. 2 is solved using the values of the parameters that 
provided the best fit to the actual temperature profile data 
points (such as those used to plot curve 1 1 0 in Figure 1 1 ). 
Figure 1 2 shows at curve 1 1 2 what could be the plot of the 
20 pure temperature perturbation of the data plotted and fitted in 
Figure 1 1 . 

[Para 44] The amplitude of the temperature effect generated 
by the injection at each valve can also be plotted, as shown in 
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Figure 1 3. In Figure 1 3, it is assumed tPiat tPiere are tPiree 
valves, a shallow valve, a deep valve, and a medium valve 
located between the shallow and deep valve. The straight line 
curve 1 1 6 of Figure 1 3 shows that no temperature effect is 
5 occurring at the shallow valve and therefore it is likely that no 
injection is occurring at such valve. The dotted line curve 1 1 8 
of Figure 1 3 shows that some temperature effect is occurring 
intermittently at the medium valve. The dashed line curve 1 20 
of Figure 1 3 shows that a temperature effect greater than that 

1 0 of the medium valve is intermittently occurring at the deep 
valve. The deep and medium valves are the more important 
valves, since the shallow valve shows no temperature effect. 
Generally, the amplitude of the perturbation is highest when 
production fluid is stationary because the gas cools the same 

1 5 fluid over a long period of time. Consequently, amplitude 
alone is a poor indicator of injection rate. 

[Para 45] The downstream decay length for each valve is also 
a useful illustration for an operator. Typically, if the 
production fluid is moving, it carries the temperature 
20 perturbation up the well. The distance that the perturbation 
persists before disappearing increases with increasing flow 
rate. Figure 1 4 shows an example plot of decay lengths versus 
time for each of the valves: shallow valve (straight line curve 
1 22) medium valve (dotted line curve 1 24), and deep valve 
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(dashed line curve 1 26). Note that the decay length increases 
at the beginning of each cycle for the two flowing valves, 
reaches a maximum, and stays at the maximum briefly before 
injection stops. 

5 [Para 46] As previously stated, the amplitude alone (see 
Figure 1 3) is a poor indication of injection rate. A better 
qualitative indicator of gas injection rate is the product of 
amplitude and downstream decay length. As with the flow 
rate, the injection rate increases early in the cycle. However, 
1 0 the injection rate starts to decline before the end of the cycle. 
Figure 1 5 illustrates a plot of the product of amplitude and 
downstream decay length versus time for each of the valves: 
shallow valve (straight line curve 1 28) medium valve (dotted 
line curve 1 30), and deep valve (dashed line curve 1 32). 

1 5 [Para 47] Contour plots can be particularly useful for an 

operator to analyze the performance of each valve at a time. 
Figure 1 6 shows such a contour plot of one of the valves, for 
instance the deep valve. With a contour plot that plots 
temperature versus time and measured depth, an operator can 

20 analyze where and when injection is occurring at a particular 
valve, as well as the extent of such injection. 

[Para 48] In operation, data from the distributed temperature 
sensor 20 is sent to the processor 22 via the communication 
link 30. The processor 20, which is loaded with the relevant 
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algorithm or model, analyzes the temperature profile data 
itself to minimize or remove any non-relevant temperature 
"noise" and to focus on the data points or sections that contain 
valuable information. Representative algorithms are illustrated 
5 in Figures 3, 5, and 6-10. The inclusion of the processor 22 
minimizes operator involvement in the analysis. Depending on 
the embodiment of the algorithm or model used in the 
processor 22, the user may be prompted by the processor 22 
to answer certain questions, such as the approximate location 
1 0 of gas lift valves. The processor 22 then presents the results 
of the analysis to the user so as to highlight the valuable 
information that was extracted from the data by the processor 
22. Examples of useful plot presentations are shown in 
Figures 11-16. 

1 5 [Para 49] The use of the present invention in relation to gas- 
lift systems, such as the one shown in Figure 4, can be 
particularly beneficial. The present invention enables an 
operator to determine the location, time, and extent of gas 
injection (i.e. the location and operation of gas lift valves) in a 

20 wellbore. 

[Para 50] Having the results of the present invention on hand, 
an operator can then diagnose problems with the well, such as 
leaking or non-operating valves or valves with sub-optimal 
characteristics. 
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[Para 51] Although the gas-lift operation was described, it is 
understood that the present invention may be used for other 
types of operations, such as identification of cross flow 
between reservoir intervals at different reservoir pressures 
5 when the well is shut in, identification of gas inflow from the 
formation through perforated intervals, wellbore 
communication investigation, steam floods, water profiles, 
optimizing sampling processes and timing, and determining 
fracture height. 

1 0 [Para 52] As previously described, instructions of the various 
routines discussed herein (such as the method and algorithm 
performed by the processor 22 and subparts thereof including 
equations and plots) may comprise software routines that are 
stored on memory 34 and loaded for execution on the CPU 32. 

1 5 Data and instructions (relating to the various routines and 
inputted data) are stored in the memory 34. The memory 34 
may include semiconductor memory devices such as dynamic 
or static random access memories (DRAMs or SRAMs), erasable 
and programmable read-only memories (EPROMs), electrically 

20 erasable and programmable read-only memories (EEPROMs) 
and flash memories; magnetic disks such as fixed, floppy and 
removable disks; other magnetic media including tape; and 
optical media such as compact disks (CDs) or digital video 
disks (DVDs). 
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[Para 53] In an alternate embodiment, the algorithm 
discussed above is modified to further identify the depths of 
injection valves or other gas-injection events rather than using 
known injection valve depths. Specifically, a match filter is 
5 constructed with the shape that is characteristic of a gas- 
injection event. The algorithm processes temperature profiles 
to identify candidate depths where gas appears to be injected. 
The algorithm discussed above (see, for example. Figure 7 and 
its associated description), can then be used to process the 
1 0 profiles at the candidate depths. However, minor 

modifications to the previously discussed algorithm can be 
made, as discussed below. 

[Para 54] The new algorithm builds a match filter having 
characteristics expected from gas injection, e.g. a sharp 
1 5 change in temperature on the upstream (deep) side of the 

event and a more gradual decay on the downstream side. The 
mathematical convolution of the filter with a profile indicates 
candidate depths. 

[Para 55] It should be noted that standard forms of matched 
20 filters may be modified to accommodate use with distributed 
temperature sensor temperature profiles. Normally, matched 
filters maximize the output signal-to-noise ratio in a filtering 
system when the noise satisfies certain characteristics, and the 
most important requirement is that the power spectrum of the 
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noise be independent of frequency. Distributed temperature 
sensor profiles can violate this requirement. In such systems, 
the dominant part of the noise tends to be the background 
trend that varies slowly In space. Consequently, the spectrum 
5 of the noise Is Inversely proportional to spatial frequency at 
low frequencies. To suppress the baclcground trend, terms can 
be added to the filter to make It orthogonal to the background. 
In one embodiment, constant and linear terms make the filter 
orthogonal to linear background trends. In this example, a 
1 0 final modification to the filter Is normalization. The amplitude 
of the convolution should be unity when the profile has an 
injection signature with unit amplitude. 

[Para 56] With the addition of the match filter, an 
identification algorithm is Illustrated in flow chart form in 

1 5 Figure 1 7. The identification algorithm is similar to the 

algorithm described above with reference to Figures 7-1 1 with 
several modifications. In an Initial step 140, temperature 
profile data Is read by the system, e.g. a system processor, 
when the data is received from, for example, a distributed 

20 temperature sensor. In a next step 142, the system processes 
convolution C of the match filter with a temperature profile. In 
a next step 144, an estimate Is made of the statistical noise 
level for each candidate depth. The model Is then fit to the 
actual data for each candidate depth, as set forth In step 146. 
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In a next step 148, the results for each candidate depth are 
tested for statistical significance. Steps 142, 144, 146 and 
148 are repeated for each profile, as set forth in step 1 50. 
Subsequently, in step 1 52, the results may be displayed in one 
5 or more ways that are beneficial and/or valuable to the user. 

[Para 57] Many of these steps have been described above with 
reference to Figures 7-1 1 , but there are several differences. 
First, processing convolution C of the match filter (step 142) is 
different. In the previous algorithm, all temperature profiles 

1 0 are processed for a particular injection depth at one time. In 
the algorithm discussed with reference to Figure 1 7, the 
candidate depths vary from temperature profile to temperature 
profile, and each profile is processed individually. The 
convolution step for a temperature profile can be described 

1 5 with reference to the flow chart illustrated in Figure 1 8. 

[Para 58] Initially, in step 1 54, the system computes 
convolution C of the match filter with a temperature profile. It 
should be noted that at shallow depths, temperature profiles 
often have significant anomalies. Accordingly, the algorithm 
20 may be designed to ignore initial distances, e.g. the first 500 
meters of depth. Although the convolution smooths the data, 
point-to-point fluctuations may still be too large. Accordingly, 
C may be further smoothed with, for example, a Savitzky- 
Colay filter (see W. H. Press et al.. Numerical Recipes in C, 2nd 



30 



Ed., page 650, Cambridge University Press, New Yorl< (1992)), 
as Illustrated in step 1 56. 

[Para 59] In a next step 1 58, local extrema are located where 
the first derivative of the smoothed convolution changes sign. 
5 When the filter is normalized as described above, the 
convolution with a cooling event is negative. Thus, local 
minima, where the second derivative is positive, are selected 
from the extrema. A threshold test is applied. For example, 
the magnitude of the convolution must exceed a threshold for 

1 0 a particular minimum to be accepted, and the convolution 
must increase by another threshold in the vicinity of the 
minimum. The number of minima that satisfy the threshold 
tests is usually small. If there are too many minima, the 
system selects minima having the largest second derivative of 

15 C, as set forth in step 1 60. Those with a smaller second 
derivative are eliminated. 

[Para 60] In a subsequent step 1 62, a mean position is used 
for multiple minima that are too close. Specifically, if minima 
occur too close to one another, the procedure for fitting 
20 multiple injection candidates simultaneously may not 

converge. Thus, when the separation of a group of candidates 
is too small, a single candidate at the mean depth replaces the 
group. The depths of the minima determine the injection 
candidates, as set forth in step 1 64. The minima that pass all 
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the tests are the injection candidates that undergo further 
processing via the algorithm illustrated in Figure 1 7. 

[Para 61 ] Referring again to Figure 1 7, subsequent to 
processing convolution C in step 142, the system estimates 
5 the statistical noise level for each candidate depth in step 1 44. 

[Para 62] Algorithm steps 146 and 148 can be performed 
similar to steps 78 and 80 described above and illustrated in 
Figure 7. Furthermore, the results may be displayed in, for 
example, graphical form and showing ranges of depths in 
1 0 which events are clustered. 

[Para 63] The present invention may be used with land as well 
as subsea wellbores, including subsea wellbores with 
subsurface gas-lift installations. 

[Para 64] Moreover, the results of the present invention may 
1 5 be combined with other measurements to analyze a well's 
performance more thoroughly and to help decide how to 
improve performance. For instance, the present invention may 
be combined with measurements of flow rate or pressure. 

[Para 65] As previously described, instructions of the various 
20 routines discussed herein (such as the method and algorithm 
performed by the processor 22 and subparts thereof including 
equations and plots) may comprise software routines that are 
stored on memory 34 and loaded for execution on the CPU 32. 
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Data and instructions (relating to tine various routines and 
inputted data) are stored in the memory 34. The memory 34 
may include semiconductor memory devices such as dynamic 
or static random access memories (DRAMs or SRAMs), erasable 
5 and programmable read-only memories (EPROMs), electrically 
erasable and programmable read-only memories (EEPROMs) 
and flash memories; magnetic disks such as fixed, floppy and 
removable disks; other magnetic media including tape; and 
optical media such as compact disks (CDs) or digital video 
10 disks (DVDs). 

[Para 66] In the following embodiment, models/algorithms 
are provided for indicating a flow rate of a fluid produced 
through production tubing 16. As illustrated in Figure 19, the 
fluid flow rate is modeled via a specific well model relating 

1 5 temperature characteristics to the flow rate through tubing 1 6 
(see block 166). Once the model is established, temperatures 
may be measured along the well (see block 1 68) with, for 
example, distributed temperature sensor 20. The measured 
temperatures are applied to the well model (see block 1 70) 

20 which utilizes the temperature characteristics to determine a 
fluid flow rate (see block 1 72). 

[Para 67] The algorithms discussed above with reference to 
Figures 1-16 can be used to identify where and when gas is 
injected into tubing 16 from the casing/tubing annulus of a 
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gas-lift well. Those algorithms fit an exponential function to 
the spatial distribution of the temperature perturbation caused 
by gas injection. As is discussed below, however, the 
downstream decay length of the exponential is related to the 
5 flow rate in tubing 1 6 and to the rate of radial heat transport 
between the tubing and the surrounding formation. 
Furthermore, the amplitude of the temperature perturbation 
also can be related to the flow rate. 

[Para 68] Gas-lift injection modifies the temperature profile 
10 along a gas-lift well. Both the amplitude and the shape of the 
perturbation depend on the production fluid flow rate. 
Accordingly, the temperature perturbation can be measured 
by, for example, distributed temperature sensor 20 and used 
to determine flow rate. An automated process of determining 
1 5 such flow rates is illustrated generally in Figures 20 and 21 . 

[Para 69] Referring first to Figure 20, flow rate can be 
determined based on the downstream decay length of the 
temperature perturbation. In this example, temperature data 
from sensor 20 is input to processor 22, as illustrated in block 
20 1 74. The temperature data enables evaluation of the decay 
length of the thermal perturbation at a gas injection location, 
as illustrated by block 1 76. Once the decay length is 
determined, a model relating decay length and flow rate is 
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applied to determine tlie flow rate of production fluid through 
tubing 1 6, as Illustrated by block 1 78. 

[Para 70] Similarly and with reference to Figure 21 , flow rate 
also can be determined based on the amplitude of the 
5 temperature perturbation. Again, temperature data from, for 
example, distributed temperature sensor 20 is input to 
processor 22, as illustrated by block 1 80. The temperature 
data enables determination of the amplitude of the thermal 
perturbation at a given a gas injection location, as illustrated 
1 0 by block 1 82. The thermal perturbation data can then be 
utilized by a model relating amplitude to flow rate to 
determine the flow rate of production fluid through tubing 1 6, 
as illustrated by block 1 84. 

[Para 71 ] Examples of specific models that can be used to 
1 5 determine the production fluid flow rates are discussed in 

detail below. However, it should be noted that the processing 
of measured thermal data according to the models may be 
carried out on processor 22 or other suitable processing 
system. Similarly, the mathematical models/algorithms can be 
20 stored, for example, at memory 34 or other suitable location. 

[Para 72] In applying the model or models to a given set of 
thermal data, other well related parameters may be 
incorporated into the modeling to improve the accuracy of the 
determined flow rates based on the temperature profile. The 
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desirability of incorporating sucii parameters into application 
of the model may depend on such factors as gas-lift well 
environment and gas-lift system design. 

[Para 73] Referring generally to Figure 22, a variety of well 
5 related parameters 1 86, 1 88, 1 90, 1 92 and 1 94 can be utilized 
to improve the accuracy of the results when applying a given 
model, as illustrated by block 1 96. Examples of such 
parameters comprise heat capacity of the production fluid 1 86, 
thermal conductivity of the surrounding formation 1 88, 

1 0 thermal history of the well 1 90, radial heat transport in the 

well in the surrounding formation 192 (particularly when using 
the model relating decay length and flow rate) and pressure 
drop 1 94 between the annulus and tubing 1 6 (particularly 
when using the model relating thermal perturbation amplitude 

1 5 and flow rate). 

[Para 74] Whether estimating the flow rate based on the decay 
length or the amplitude of the thermal perturbation, produced 
fluid heat capacity is a parameter that often affects 
quantitative estimates of the production fluid flow rate. It 
20 should be noted, however, that the produced fluid can be a 
mixture of fluids, such as water and oil. The heat capacity per 
unit mass of water is typically three times as large as that of 
oil. Consequently, uncertainty in the produced-water fraction 
causes an equal or larger relative uncertainty in an estimated 
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flow rate. If the produced-water fraction is known from 
surface measurements, the heat capacity of the fluid may be 
determined or estimated. Otherwise, the water fraction of the 
produced fluid can be measured. In some applications, a 
5 differential pressure measurement immediately below the gas- 
injection depth can be used to provide the water fraction, 
assuming there is no gas at that point. 

[Para 75] In the following discussion, embodiments of models 
are discussed and developed to facilitate an understanding of 

1 0 the ability to determine flow rates based on temperature 

profiles in gas-lift wells, as graphically illustrated in Figures 20 
and 21 . In these examples, certain assumptions are made 
about the design of the gas-lift well. Specifically, a vertical 
well is considered with fluid entering production tubing 16 at 

1 5 the local geothermal temperature. Injected gas flows down the 
casing/tubing annulus and enters the bottom of tubing 16. 
Production fluid and the injected gas flow to a surface location 
through production tubing 16. Furthermore, it is assumed that 
the axial and radial heat transport in the tubing and annulus 

20 are in steady-state, but the radial heat transport in the 
formation is time dependent. The model can be modified 
easily to account for an inclined well or gas injection above the 
bottom of the tubing. 
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[Para 76] The models utilized are based on heat transport 
equations and their solutions. A solution strategy is to 
compute the net axial transport of enthalpy into small 
segments of the production tubing and annulus and to equate 
5 these to the net losses from the segments by radial heat 
transport. The mathematical basis is established as follows: 

Equation 3 



1 0 The left side of the equation is axial transport of enthalpy, the 
right side is net radial heat loss. If the continuous Joule- 
Thompson OT) effect is ignored, the enthalpy may be replaced 
with the heat capacity: 
Equation 4 



The heat flow rate from the annulus to the formation in the 

interval dz is: 
Equation 5 



dz dz 



15 




1 d(qp^-qj 
dz 



<lFa=^aC,a(rG-TJdzlA, 



20 




38 



The dimensionless time, r , accounts for the difference 
between the local geothermal temperature, Te , and the actual 
temperature of the formation at the well. The heat flow rate 
from the tubing to the annulus in the interval dz is: 
5 Equation 6 

qa,=w^c^^(J^-T,)dzlB, 



Combining Eq. 4-6, the following is obtained: 
Equation 7 

10 — — — ^ 

dz A B 
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In the tubing, a similar derivation results in: 
Equation 8 



dT, _ T,-T^ 
~dl~ B' 



Elimination of 71 in Eqs. 7 and 8 produces a second order 
differential equation for the tubing temperature: 
Equation 9 

AB'^-^ + B"^-T+ T = 0, 

dz' dz ' " 
B'' = B' + AB'/B-A 
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For a linear geothermal gradient, the solution to Eq. 9 is: 
Equation 10 

= (1 - \B')ae''' + (1 -X,B')Pe''' +g^{B' + z- B') + r„ , 
\ = (ylB'^+4AB' - B') 1 2AB\ 
= -{B' + Vfi'^ + AAB') 1 2AB'. 

5 

It should be noted that Ai is positive, and terms involving it 
are usually important only at the bottom of the well. X 2 is 
negative, and terms involving it are usually important only near 
the surface. The boundary conditions are the inlet 
1 0 temperature of the gas at the surface and the tubing 
temperature at bottom hole: 
Equation 1 1 

T^(z = 0)=T^, 



1 5 The cooling coefficient, R, adds the JT effect at the gas- 
injection point. Evaluating and T^Mfrom Eq. 10, it can be 
determined: 
Equation 12 
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(l-A.B'XGggB'+DG') 
(l-Z,B')G'-(l-AiB')G'' 
(l-A^B')DG' + (l-A,B')GgaB' 
(l-A^B')G'-(l-Z^B')G' 
j^^ T^-gG(B'-B')-T^, 
l-A,B' 

^^ w^iR+c^^g^B') + w^c^g^B' 



eye can be neglected. In such approximation, the 
constants simplify to: 
Equation 1 3 



GggB' ^ w^iR + c^^g^B') + w^c^g^B'' 



(l-Ji^B )DG' + (1 - X^B')Gg(.B' 
{\-X,B)G' 



[Para 77] Heat transfer coefficients can be estimated from the 
dimenslonless Nusselt number: 




Nud =Udlk 
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In laminar flow conditions in a pipe, Nud is 4.4, and in single- 
phase turbulent flow, Nud may be estimated using the 
Reynolds number and the Prandtl number as follows: 
Equation 14 

Nud ^O.OISRCd'''^ Pr" 
Pr^cju/k 

When the pipe is warmer than the fluid, n is 0.3, otherwise n is 
0.4. Also, in the annulus, the diameter cf '\s replaced by the 
hydraulic diameter 2(rc-rt). 

1 0 [Para 78] By way of example, methane can be used as an 

injection fluid with gas-lift system 1 0. Because the viscosity of 
methane is small, the Reynolds number is usually greater than 
1 0,000. Also, because multi-phase flow enhances turbulence, 
radial transport in the production tubing is expected to be very 

1 5 efficient. The tubing heat-transfer coefficient is assumed to be 
much greater than the annulus coefficient. Consequently, the 
heat transfer between the tubing and the annulus involves only 
annulus properties. Furthermore, the coefficient for heat 
transfer between the tubing and the annulus is assumed to 

20 equal the coefficient for transfer between the annulus and the 
formation. That is: 

Equation 1 5 
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[Para 79] In some applications, the mathematical basis of the 
models can be simplified. Consider first the ratio C'/C". It is 
5 proportional to: 

The other factors in the ratio are of order unity, and the 
exponent is: 

-Jb'^+aab' liab' 

1 0 The exponent is smallest when ivp can be neglected and iva is 
large. In that case, the first term in the numerator of >4 can be 
neglected and the exponent of C'/C". becomes: 

Vl/A'+4/A5 L 

1 5 C"can be neglected when the exponent is greater than 2%, i.e., 
when: 

Equation 16 

20 In typical cases, C'can be neglected when the gas flow rate of, 
for example, methane is less than 5 kg/s. 
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[Para 80] For clarification, tlie nomenclature used herein is as 
follows: 



Symbol 


Meaning 


Units 
(SI) 




Symbol 


Meaning 


Units 
(SI) 


c 


Heat 

capacity 


j/(kg 
K) 






Mass flow rate 


kg/s 


d 


Diameter 


m 






Depth 


m 


Or 


Geothermal 
gradient 


K/m 










H 


Enthalpy per 
unit mass 


i/ka 




X 


Decay rate 


1 /m 


k 


Thermal 
conductivity 


w/(Km) 






Viscosity 


Pas 


L 


Gas- 
injection 
depth 


m 




V 


Mean velocity 


m/s 


q 


Heat 

transfer rate 


w 




P 


Density 


kg/m3 


r 


Radius 


m 










R 


JT cooling 
coefficient 


j/kg 




Nud 


Nusselt 
number 


none 


t 


Time 


s 




Pr 


Prandtl 
number 


none 
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T 


Temperature 


K 






Reynolds 
number 


none 


U 


Heat 


w/K 




T 


Dimensionless 


none 




transfer 








time 






coefficient 













[Para 81 ] Furtliermore, tlie definition of tlie various subscripts 
is as follows: 

a-annulus; M-bottom hole; c-casing; e-Earth; F-formation; g- 
5 gas; C"-geothermal; p-production fluid; 5-surface of Earth; t- 
inside tubing. 

[Para 82] To estimate flow rate from the amplitude, the 
thermal discontinuity is first determined from the temperature 
profile and then the following equation is solved for flow rate: 

1 0 Equation 1 7 

In this example, the second term in Eq. 10 for 7"f has been 
1 5 neglected, because the exponential factor suppresses it near 
the bottom of the well. In many cases, the second terms in the 
numerator and the denominator of Eq. 1 7 are much smaller 
than the first terms. The discontinuity is approximately equal 
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to the total cooling power divided by the flow rate and heat 
capacity of the production fluid. The effect of the heat 
capacity of the gas is reduced in Eq. 1 7 because gas is cooled 
as it approaches the injection point from above. It should 
5 further be noted the solution of Eq. 1 7 is possible when all 
injected gas is injected through a single orifice and the total 
gas flow rate is known. In this application, the solution is 
insensitive to Earth properties and thermal history. However, 
the cooling coefficient, which depends on the gas properties 

1 0 and the pressure difference between the annulus and the 

tubing at the injection depth, is needed for the solution. The 
pressure change in the annulus from the surface to the 
injection depth is small, because the gas density is 
comparatively low and the frictional pressure gradient 

1 5 counteracts the gravitational gradient. Thus, the annulus 

pressure at depth may be estimated accurately, and the tubing 
pressure is measured. 

[Para 83] To determine the flow rate from the decay length of 
the thermal perturbation produced by gas injection, the decay 
20 length is first obtained from the temperature profile. Then, 
flow rate may be determined by solving for X i of Eq. 1 0. 
However, the solution depends on several parameters, 
including heat capacity of the produced fluid and radial heat 
transport in the well. Additionally, the dimensionless time r 
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in the parameter A depends on the earth's thermal diffusivity 
and the thermal history of the well. An approximation to the 
analytical solution of the diffusion equation uses a constant 
heat flux. When the time t is much greater than p^,rc/ke 
5 (typically a few hours), analytical solutions for different 
boundary conditions become indistinguishable. Therefore, 
details of distant thermal history of the well can be ignored, 
but recent history can be important. Accurate flow-rate 
estimates can benefit from a numerical solution of the 
1 0 diffusion equation with the measured temperature history as 
the boundary condition. 

[Para 84] By way of further explanation, the decay length of 
the Joule-Thompson cooling perturbation is 1 / Ai. During 
normal production, the total heat capacity of the production 
1 5 fluid is much greater than the total heat capacity of the 
injected gas, i.e., B'»B. In this approximation, the decay 
length is: 

Equation 1 8 



20 



The first term is the effect of heat transfer between the tubing 
and the annulus. The other terms are the effect of heat 
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transfer between the annulus and the formation. Important 
factors are the total heat capacity of the production fluid, the 
heat transfer coefficients and the dimensionless time. 

[Para 85] Accordingly, the mathematical models discussed 
5 above can be used to determine flow rates in a gas-lift well 
based on either or both the decay length and the amplitude of 
the injection-induced thermal perturbation. Other parameters 
also can be useful in improving the accuracy of the determined 
flow rate. For example, heat capacity of the production fluid 

1 0 can be important when relying on either decay length models 
or amplitude models. When using an amplitude model, it can 
be important to measure the pressure drop between the 
annulus and the tubing. When using a decay length model, it 
often is helpful to determine the radial heat transport in both 

1 5 the well and the surrounding formation. Furthermore, in the 
embodiments described, the temperature data collection, 
application of a model to the temperature data, and the 
determination of flow rates are conducted on processor system 
22. A variety of a graphical displays or other output formats 

20 may be displayed on output device 38 to convey flow rate 
information to a system operator. 

[Para 86] In another embodiment, the gas lift performance of 
a well can be optimized by utilizing the downhole fluid flow 
rates determined through the temperature data obtained, for 
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example, via distributed temperature sensor 20. Gas injected 
into many gas-lift applications is not within an optimal range 
due to, for example, operators injecting too much gas into the 
wellbore. The algorithms discussed above for determining 
5 flow rate in a gas-injection well can be used in the present 
embodiment to determine fluid flow rates. Additionally, the 
algorithms can be adjusted to provide a feedback loop that 
enables automatic changes to the gas injection rate and 
computation of the optimal amount of gas injection to 
1 0 maximize the fluid flow rate. 

[Para 87] Referring generally to Figure 23, the general 
methodology is illustrated in flow chart form. Initially, a flow 
rate of the produced fluid is determined based on one or more 
wellbore parameters (see block 198). Subsequently, an 

1 5 analysis is performed as to whether the flow rate is in an 
optimal range (see block 200). In this embodiment, the 
analysis is performed automatically via, for example, processor 
22. The optimal range can be determined in a variety of ways, 
including use of data from similar wells, use of historical data 

20 from the well being analyzed or by adjusting the gas injection 
rate and tracking whether the production fluid flow rate is 
increasing or decreasing. If the processor determines the flow 
rate is not optimized, an action is taken, e.g. changing the gas 
injection rate, to adjust the flow rate (see block 202). 
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Following adjustment, the new fluid flow rate is again 
determined (see block 1 98) and the process Is repeated. 

[Para 88] As discussed above with reference to Figures 1 9-22, 
the flow rate of fluid produced through production tubing 16 
5 can be obtained from temperatures measured along the 
wellbore. As illustrated in Figure 24, the fluid flow rate Is 
modeled via a specific well model/algorithm that relates 
temperature to the production fluid flow rate through tubing 
1 6 (see block 204). After establishing the suitable model, 

1 0 temperatures are measured along the well (see block 206). A 
distributed temperature sensor, such as the distributed 
temperature sensor 20 discussed above, works well to obtain a 
temperature profile that can automatically be provided to 
processor 22. The measured temperatures are applied to the 

1 5 well model (see block 208) which uses those measured 

temperatures to determine a fluid flow rate of the production 
fluid (see block 210). In this embodiment, however, the 
model/algorithm is expanded to automatically optimize that 
fluid flow rate. 

20 [Para 89] Processor 22 can be used in a closed loop feedback 
system to facilitate this flow rate optimization by continually 
analyzing whether the flow rate is within a determined optimal 
range. Specifically, upon determining a fluid flow rate, the 
algorithm performs a first test 21 2 and checks to see if the 
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fluid flow rate through the production tubing is too fast, e.g. 
above the optimal range. If the flow rate is too fast, processor 
22 acts to decrease the fluid flow rate (see block 214) by, for 
example, decreasing the flow of injection gas. The process is 
5 then returned to block 206 to once again measure 

temperatures along the well for determining the new flow rate. 
If, however, the first test 212 does not detect a fluid flow that 
is too fast, a second test 216 checks to see if the flow rate is 
too slow. If the flow rate is too slow, processor 22 acts to 

1 0 increase the fluid flow rate (see block 21 8) by, for example, 
increasing the gas injected. The process is then returned to 
block 206 to again measure temperatures along the well for 
determining the new flow rate. When second test 21 6 is 
performed and the fluid flow rate in the production tubing is 

1 5 not too slow, then the flow rate is in the optimal range and the 
process is returned to block 206 for subsequent checking of 
the fluid flow rate. Thus, use of the algorithms discussed 
above can be automated to continually check and optimize the 
production fluid flow rate. 

20 [Para 90] While the present invention has been described with 
respect to a limited number of embodiments, those skilled in 
the art, having the benefit of this disclosure, will appreciate 
numerous modifications and variations therefrom. It is 
intended that the appended claims cover all such modifications 
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and variations as fall within the true spirit and scope of this 
present Invention. 
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